Phylogenomics of the world's otters.

Comparative whole-genome analyses hold great power to illuminate commonalities and differences in the evolution of related species that share similar ecologies. The mustelid subfamily Lutrinae includes 13 currently recognized extant species of otters,1-5 a semiaquatic group whose evolutionary history is incompletely understood. We assembled a dataset comprising 24 genomes from all living otter species, 14 of which were newly sequenced. We used this dataset to infer phylogenetic relationships and divergence times, to characterize patterns of genome-wide genealogical discordance, and to investigate demographic history and current genomic diversity. We found that genera Lutra, Aonyx, Amblonyx, and Lutrogale form a coherent clade that should be synonymized under Lutra, simplifying the taxonomic structure of the subfamily. The poorly known tropical African Aonyx congicus and the more widespread Aonyx capensis were found to be reciprocally monophyletic (having diverged 440,000 years ago), supporting the validity of the former as a distinct species. We observed variable changes in effective population sizes over time among otters within and among continents, although several species showed similar trends of expansions and declines during the last 100,000 years. This has led to different levels of genomic diversity assessed by overall heterozygosity, genome-wide SNV density, and run of homozygosity burden. Interestingly, there were cases in which diversity metrics were consistent with the current threat status (mostly based on census size), highlighting the potential of genomic data for conservation assessment. Overall, our results shed light on otter evolutionary history and provide a framework for further in-depth comparative genomic studies targeting this group.

The role and impact of Zootaxa in mammalogy in its first 20 years.

Zootaxa came as a new and innovative publication medium for taxonomy, amidst a scenario of devaluation of this important biological science. After 20 years, it has ascertained itself as one of the main journals in animal taxonomy. However, the contribution of the journal to the taxonomy of Mammalia (mammals), one of the most studied groups of animals with a long-standing, dedicated spectrum of specialized journals (mammalogy), could have been expected as minor. All the current and former editors of the Mammalia section of Zootaxa analyzed the relative contribution of the journal to the description of new species of mammals since 2001. We also analyzed the contribution of Zootaxa by taxon, geographic origin of taxa, and geographic origin of first authors. The taxonomic methodology of authors in species description is described as well as the temporal trends in publications and publication subjects. We highlight the editors' picks and eventually, the challenges for the future. We found that Zootaxa has had a significant contribution to mammalogy, being the second journal (the first being Journal of Mammalogy) in terms of number of new species described (76; 10.6% of the new mammalian species described between 2001 and 2020). The majority of the new species were described following an integrative taxonomic approach with at least two sources of data (86%). The analysis of published taxa, their geographic origin, and the country of origin of first authors shows a wide coverage and exhaustive representation, except for the species from the Nearctic. We conclude that Zootaxa has likely responded to a repressed demand for an additional taxonomic journal in mammalogy, with as possible appeals the absence of publication fees and an established publication speed. With 246 articles published in the past 20 years, the Mammalia section of Zootaxa embraces a large spectrum of systematic subjects going beyond alpha taxonomy. The challenges for the future are to encourage publications of authors from the African continent, still poorly represented, and from the palaeontology community, as the journal has been open to palaeontology since its early days.

Correction: Evolutionary history of Carnivora (Mammalia, Laurasiatheria) inferred from mitochondrial genomes.

[This corrects the article DOI: 10.1371/journal.pone.0240770.].

Evolutionary history of Carnivora (Mammalia, Laurasiatheria) inferred from mitochondrial genomes.

The order Carnivora, which currently includes 296 species classified into 16 families, is distributed across all continents. The phylogeny and the timing of diversification of members of the order are still a matter of debate. Here, complete mitochondrial genomes were analysed to reconstruct the phylogenetic relationships and to estimate divergence times among species of Carnivora. We assembled 51 new mitogenomes from 13 families, and aligned them with available mitogenomes by selecting only those showing more than 1% of nucleotide divergence and excluding those suspected to be of low-quality or from misidentified taxa. Our final alignment included 220 taxa representing 2,442 mitogenomes. Our analyses led to a robust resolution of suprafamilial and intrafamilial relationships. We identified 21 fossil calibration points to estimate a molecular timescale for carnivorans. According to our divergence time estimates, crown carnivorans appeared during or just after the Early Eocene Climatic Optimum; all major groups of Caniformia (Cynoidea/Arctoidea; Ursidae; Musteloidea/Pinnipedia) diverged from each other during the Eocene, while all major groups of Feliformia (Nandiniidae; Feloidea; Viverroidea) diversified more recently during the Oligocene, with a basal divergence of Nandinia at the Eocene/Oligocene transition; intrafamilial divergences occurred during the Miocene, except for the Procyonidae, as Potos separated from other genera during the Oligocene.

Phenotypic integration in feliform carnivores: Covariation patterns and disparity in hypercarnivores versus generalists.

The skeleton is a complex arrangement of anatomical structures that covary to various degrees depending on both intrinsic and extrinsic factors. Among the Feliformia, many species are characterized by predator lifestyles providing a unique opportunity to investigate the impact of highly specialized hypercarnivorous diet on phenotypic integration and shape diversity. To do so, we compared the shape of the skull, mandible, humerus, and femur of species in relation to their feeding strategies (hypercarnivorous vs. generalist species) and prey preference (predators of small vs. large prey) using three-dimensional geometric morphometric techniques. Our results highlight different degrees of morphological integration in the Feliformia depending on the functional implication of the anatomical structure, with an overall higher covariation of structures in hypercarnivorous species. The skull and the forelimb are not integrated in generalist species, whereas they are integrated in hypercarnivores. These results can potentially be explained by the different feeding strategies of these species. Contrary to our expectations, hypercarnivores display a higher disparity for the skull than generalist species. This is probably due to the fact that a specialization toward high-meat diet could be achieved through various phenotypes. Finally, humeri and femora display shape variations depending on relative prey size preference. Large species feeding on large prey tend to have robust long bones due to higher biomechanical constraints.

Author Correction: The globally invasive small Indian mongoose Urva auropunctata is likely to spread with climate change.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The globally invasive small Indian mongoose Urva auropunctata is likely to spread with climate change.

Invasive alien species represent one of the major factors of global loss of biodiversity and disruption of natural ecosystems. The small Indian mongoose, Urva auropunctata, is considered one of the wild carnivore species with the greatest negative impact on global biodiversity. Understanding of the factors underpinning the species' distribution and potential dispersion in a context of climate change thus appears crucial in the conservation of native ecosystems. Here we modelled the current and future climatically favourable areas for the small Indian mongoose using Ecological Niche Modelling based on data sets filtrated in environmental spaces. Projections from these models show extensive current favourable geographical areas, covering continental and insular regions within tropical and sub-tropical latitudes. Moreover, predictions for 2050 reveal that climate change is likely to expand current favourable areas north of the current favourable spaces, particularly in Eastern Europe. This climate-induced expansion is particularly worrisome given that the species is already spreading in the Balkan region. Our projections suggest that it is very likely that the small Indian mongoose will have an increasing influence on ecosystems and biodiversity in Europe by 2050.

Current and future climatic regions favourable for a globally introduced wild carnivore, the raccoon Procyon lotor.

Invasive species are considered as one of the major threats to biodiversity and represent a major challenge in the conservation of natural ecosystems, in preventing damage to agricultural production, and human health risks. Environmental Niche Modelling has emerged as a powerful tool to predict the patterns of range expansion of non-native species and to direct effective strategies for managing biological invasions. The raccoon, Procyon lotor, is a wild mesocarnivore presenting a high adaptability and showing successful introduced populations worldwide. Here, we modelled the current and future climatically favourable areas for the raccoon using two protocols, based on data sets filtrated in geographic and environmental spaces. Projections from these models show extensive current favourable geographical areas covering extensive regions of temperate biomes. Moreover, predictions for 2050 reveals extensive new favourable areas north of the current favourable regions. However, the results of the two modeling approaches differ in the extent of predicted favourable spaces. Protocols using geographically filtered data present more conservative forecasts, while protocol using environmental filtration presents forecasts across greater areas. Given the biological characteristics and the ecological requirements of a generalist carnivore such as the raccoon, the latter forecasts appears more relevant and should be privileged in the development of conservation plans for ecosystems.

Using metabarcoding to compare the suitability of two blood-feeding leech species for sampling mammalian diversity in North Borneo.

The application of high-throughput sequencing (HTS) for metabarcoding of mixed samples offers new opportunities in conservation biology. Recently, the successful detection of prey DNA from the guts of leeches has raised the possibility that these, and other blood-feeding invertebrates, might serve as useful samplers of mammals. Yet little is known about whether sympatric leech species differ in their feeding preferences, and whether this has a bearing on their relative suitability for monitoring local mammalian diversity. To address these questions, we collected spatially matched samples of two congeneric leech species Haemadipsa picta and Haemadipsa sumatrana from lowland rainforest in Borneo. For each species, we pooled ~500 leeches into batches of 10 individuals, performed PCR to target a section of the mammalian 16S rRNA locus and undertook sequencing of amplicon libraries using an Illumina MiSeq. In total, we identified sequences from 14 mammalian genera, spanning nine families and five orders. We found greater numbers of detections, and higher diversity of OTUs, in H. picta compared with H. sumatrana, with rodents only present in the former leech species. However, comparison of samples from across the landscape revealed no significant difference in mammal community composition between the leech species. We therefore suggest that H. picta is the more suitable iDNA sampler in this degraded Bornean forest. We conclude that the choice of invertebrate sampler can influence the detectability of different mammal groups and that this should be accounted for when designing iDNA studies.

Genetic Structure and Phylogeography of the Leopard Cat (Prionailurus bengalensis) Inferred from Mitochondrial Genomes.

The Leopard cat Prionailurus bengalensis is a habitat generalist that is widely distributed across Southeast Asia. Based on morphological traits, this species has been subdivided into 12 subspecies. Thus far, there have been few molecular studies investigating intraspecific variation, and those had been limited in geographic scope. For this reason, we aimed to study the genetic structure and evolutionary history of this species across its very large distribution range in Asia. We employed both PCR-based (short mtDNA fragments, 94 samples) and high throughput sequencing based methods (whole mitochondrial genomes, 52 samples) on archival, noninvasively collected and fresh samples to investigate the distribution of intraspecific genetic variation. Our comprehensive sampling coupled with the improved resolution of a mitochondrial genome analyses provided strong support for a deep split between Mainland and Sundaic Leopard cats. Although we identified multiple haplogroups within the species' distribution, we found no matrilineal evidence for the distinction of 12 subspecies. In the context of Leopard cat biogeography, we cautiously recommend a revision of the Prionailurus bengalensis subspecific taxonomy: namely, a reduction to 4 subspecies (2 mainland and 2 Sundaic forms).

Phylogeography of the Small Indian Civet and Origin of Introductions to Western Indian Ocean Islands.

The biogeographic dynamics affecting the Indian subcontinent, East and Southeast Asia during the Plio-Pleistocene has generated complex biodiversity patterns. We assessed the molecular biogeography of the small Indian civet (Viverricula indica) through mitogenome and cytochrome b + control region sequencing of 89 historical and modern samples to (1) establish a time-calibrated phylogeography across the species' native range and (2) test introduction scenarios to western Indian Ocean islands. Bayesian phylogenetic analyses identified 3 geographic lineages (East Asia, sister-group to Southeast Asia and the Indian subcontinent + northern Indochina) diverging 3.2-2.3 million years ago (Mya), with no clear signature of past demographic expansion. Within Southeast Asia, Balinese populations separated from the rest 2.6-1.3 Mya. Western Indian Ocean populations were assigned to the Indian subcontinent + northern Indochina lineage and had the lowest mitochondrial diversity. Approximate Bayesian computation did not distinguish between single versus multiple introduction scenarios. The early diversification of the small Indian civet was likely shaped by humid periods in the Late Pliocene-Early Pleistocene that created evergreen rainforest barriers, generating areas of intra-specific endemism in the Indian subcontinent, East, and Southeast Asia. Later, Pleistocene dispersals through drier conditions in South and Southeast Asia were likely, giving rise to the species' current natural distribution. Our molecular data supported the delineation of only 4 subspecies in V. indica, including an endemic Balinese lineage. Our study also highlighted the influence of prefirst millennium AD introductions to western Indian Ocean islands, with Indian and/or Arab traders probably introducing the species for its civet oil.

The complete mitochondrial genome of the Spotted Linsang, Prionodon pardicolor, the first representative from the family Prionodontidae (Mammalia, Carnivora).

The complete mitochondrial genome of the Spotted Linsang, Prionodon pardicolor, was sequenced using overlapping PCR products. The genome is 16,718 base pairs in length and contains the 37 genes found in a typical mammalian genome: 13 protein-coding genes, 22 transfer RNA genes and 2 ribosomal RNA genes. The overall base composition on the L-strand is A: 32.4%, C: 25.0%, G: 13.9%, T: 28.7%. The control region of P. pardicolor includes both RS2 and RS3 tandem repeats. Phylogenetic analyses support a sister relationship with the Felidae.

The complete mitochondrial genome of the boky-boky, Mungotictis decemlineata, the first representative of the Malagasy carnivores (Mammalia, Carnivora, Eupleridae).

The complete mitochondrial genome of the boky-boky, Mungotictis decemlineata, was sequenced using overlapping PCRs. The genome is 16,910 base pairs in length and contains the 37 genes found in a typical mammalian genome: 13 protein-coding genes, 22 transfer RNA genes and 2 ribosomal RNA genes. The overall base composition on the L-strand is A: 32.1%, C: 27.8%, G: 14.5%, T: 25.6%. The control region of M. decemlineata includes both RS2 and RS3 tandem repeats.

Ancestral capture of syncytin-Car1, a fusogenic endogenous retroviral envelope gene involved in placentation and conserved in Carnivora.

Syncytins are envelope protein genes of retroviral origin that have been captured for a function in placentation. Two such genes have already been identified in simians, two distinct, unrelated genes have been identified in Muridae, and a fifth gene has been identified in the rabbit. Here, we searched for similar genes in the Laurasiatheria clade, which diverged from Euarchontoglires--primates, rodents, and lagomorphs--shortly after mammalian radiation (100 Mya). In silico search for envelope protein genes with full-coding capacity within the dog and cat genomes identified several candidate genes, with one common to both species that displayed placenta-specific expression, which was revealed by RT-PCR analysis of a large panel of tissues. This gene belongs to a degenerate endogenous retroviral element, with precise proviral integration at a site common to dog and cat. Cloning of the gene for an ex vivo pseudotype assay showed fusogenicity on both dog and cat cells. In situ hybridization on placenta sections from both species showed specific expression at the level of the invasive fetal villi within the placental junctional zone, where trophoblast cells fuse into a syncytiotrophoblast layer to form the maternofetal interface. Finally, we show that the gene is conserved among a series of 26 Carnivora representatives, with evidence for purifying selection and conservation of fusogenic activity. The gene is not found in the Pholidota order and, therefore, it was captured before Carnivora radiation, between 60 and 85 Mya. This gene is the oldest syncytin gene identified to date, and it is the first in a new major clade of eutherian mammals.

Molecular phylogeny of the Herpestidae (Mammalia, Carnivora) with a special emphasis on the Asian Herpestes.

Until now, phylogenetic studies of the mongooses (Carnivora, Herpestidae) have not included an exhaustive sampling of the Asian members of this family. In this study, we used mitochondrial (Cytochrome b and ND2), nuclear (beta-fibrinogen intron 7 and Transthyretin intron 1) sequences from almost all of the recognized mongoose species to produce a well-resolved phylogeny of the Herpestidae. We also performed molecular dating analyses to infer divergence dates of the different lineages within the Herpestidae. Our results confirmed the paraphyly of the Herpestes genus and other phylogenetic relationships, which previously had only been moderately supported. The Asian herpestid species were found to form a monophyletic group within the Herpestidae. Within the Asian species, a cyto-nuclear conflict was discovered between the small Indian mongoose (Herpestes auropunctatus), the Indian gray mongoose (Herpestes edwardsii) and the Javan mongoose (Herpestes javanicus), which may have occurred through interspecific hybridization. This study inferred an Early Miocene origin for the Herpestidae and a Middle Miocene origin for the Asian mongooses.

Phylogenetic relationships of the Asian palm civets (Hemigalinae & Paradoxurinae, Viverridae, Carnivora).

The Viverridae (Mammalia, Carnivora), one of the least studied groups of carnivorans, include two subfamilies of Asian palm civets: Hemigalinae and Paradoxurinae. The relationships between and within these two subfamilies have never been thoroughly tested using an extensive molecular sample set. In this study, we gathered sequences of four genes (two mitochondrial: Cytochrome b and ND2 and two nuclear: beta-fibrinogen intron 7 and IRBP exon 1) for eight of the eleven extant species representing these two subfamilies. The results showed that: (1) the Asian palm civets (Hemigalinae and Paradoxurinae) have a single origin and form the sister-group of the (Genettinae+Viverrinae) clade, (2) the Hemigalinae (including the otter civet Cynogale bennettii) are monophyletic, (3) the Paradoxurinae are monophyletic and (4) the small-toothed palm civet (Arctogalidia trivirgata) is an early offshoot within the Paradoxurinae. Using a relaxed molecular clock analysis, the differentiation of the (Hemigalinae+Paradoxurinae) was inferred to occur in the Late Oligocene/Early Miocene.

Multigene phylogeny of the Mustelidae: resolving relationships, tempo and biogeographic history of a mammalian adaptive radiation.

BACKGROUND: Adaptive radiation, the evolution of ecological and phenotypic diversity from a common ancestor, is a central concept in evolutionary biology and characterizes the evolutionary histories of many groups of organisms. One such group is the Mustelidae, the most species-rich family within the mammalian order Carnivora, encompassing 59 species classified into 22 genera. Extant mustelids display extensive ecomorphological diversity, with different lineages having evolved into an array of adaptive zones, from fossorial badgers to semi-aquatic otters. Mustelids are also widely distributed, with multiple genera found on different continents. As with other groups that have undergone adaptive radiation, resolving the phylogenetic history of mustelids presents a number of challenges because ecomorphological convergence may potentially confound morphologically based phylogenetic inferences, and because adaptive radiations often include one or more periods of rapid cladogenesis that require a large amount of data to resolve. RESULTS: We constructed a nearly complete generic-level phylogeny of the Mustelidae using a data matrix comprising 22 gene segments (approximately 12,000 base pairs) analyzed with maximum parsimony, maximum likelihood and Bayesian inference methods. We show that mustelids are consistently resolved with high nodal support into four major clades and three monotypic lineages. Using Bayesian dating techniques, we provide evidence that mustelids underwent two bursts of diversification that coincide with major paleoenvironmental and biotic changes that occurred during the Neogene and correspond with similar bursts of cladogenesis in other vertebrate groups. Biogeographical analyses indicate that most of the extant diversity of mustelids originated in Eurasia and mustelids have colonized Africa, North America and South America on multiple occasions. CONCLUSION: Combined with information from the fossil record, our phylogenetic and dating analyses suggest that mustelid diversification may have been spurred by a combination of faunal turnover events and diversification at lower trophic levels, ultimately caused by climatically driven environmental changes. Our biogeographic analyses show Eurasia as the center of origin of mustelid diversity and that mustelids in Africa, North America and South America have been assembled over time largely via dispersal, which has important implications for understanding the ecology of mustelid communities.

Mosaics of convergences and noise in morphological phylogenies: what's in a viverrid-like carnivoran?

Adaptive convergence in morphological characters has not been thoroughly investigated, and the processes by which phylogenetic relationships may be misled by morphological convergence remains unclear. We undertook a case study on the morphological evolution of viverrid-like feliformians (Nandinia, Cryptoprocta, Fossa, Eupleres, Prionodon) and built the largest morphological matrix concerning the suborder Feliformia to date. A total of 349 characters grouped into four anatomical partitions were used for all species of Viverridae and viverrid-like taxa plus representatives of the Felidae, Hyaenidae, Herpestidae, and one Malagasy mongoose. Recent molecular phylogenetic analyses suggest that viverrid-like morphotypes appeared independently at least three times during feliformian evolution. We thus used a synthetic molecular tree to assess morphological evolutionary patterns characterizing the viverrid-like taxa. We examined phylogenetic signal, convergence and noise in morphological characters using (a) tree-length distribution (g1), (b) partitioned Bremer support, (c) RI values and their distribution, (d) respective contributions of diagnostic synapomorphies at the nodes for each partition, (e) patterns of shared convergences among viverrid-like taxa and other feliformian lineages, (f) tree-length differences among alternative hypotheses, and (g) the successive removal of convergent character states from the original matrix. In addition, the lability of complex morphological structures was assessed by mapping them onto the synthetic molecular tree. The unconstrained morphological analysis yielded phylogenetic groupings that closely reflected traditional classification. The use of a synthetic molecular tree (constraint) combined with our thorough morphological investigations revealed the mosaics of convergences likely to have contributed to part of the historical uncertainty over viverrid classification. It also showed that complex morphological structures could be subjected to reversible evolutionary trends. The morphological matrix proved useful in characterizing several feliformian clades with diagnostic synapomorphies. These results support the removal from the traditionally held Viverridae of several viverrid-like taxa into three distinct families: Nandiniidae (Nandinia), Prionodontidae (Prionodon), and the newly defined Eupleridae (including Cryptoprocta, Fossa, Eupleres plus all "mongoose-like" Malagasy taxa). No clearly "phylogenetically misleading" data subsets could be identified, and the great majority of morphological convergences appeared to be nonadaptive. The multiple approaches used in this study revealed that the most disruptive element with regards to morphological phylogenetic reconstruction was noise, which blured the expression of phylogenetic signal. This study demonstrates the crucial need to consider independent (molecular) phylogenies in order to produce reliable evolutionary hypotheses and should promote a new approach to the definition of morphological characters in mammals. [Constrained analysis; convergence; evolutionary scenario; Feliformia; morphology; noise; phylogenetic signal; phylogeny; Viverridae.].

Molecular systematics and origin of sociality in mongooses (Herpestidae, Carnivora).

The Herpestidae are small terrestrial carnivores comprising 18 African and Asian genera, currently split into two subfamilies, the Herpestinae and the Galidiinae. The aim of this work was to resolve intra-familial relationships and to test the origin of sociality in the group. For this purpose we analysed sequences of the complete cytochrome b gene for 18 species of Herpestidae. The results showed that the mongooses were split into three clades: (1) the Malagasy taxa (Galidiinae and Cryptoprocta), (2) the true social mongooses and (3) the solitary mongooses, each group being also supported by morphological and chromosomal data. Our results suggested unexpected phylogenetic relationships: (1) the genus Cynictis is included in the solitary mongoose clade, (2) the genera Liberiictis and Mungos are sister-group, and (3) the genus Herpestes is polyphyletic. We examined the evolution of the sociality in mongooses by combining behavioural traits with the cytochrome b data. Some of the behavioural traits provided good synapomorphies for characterizing the social species clade, showing the potential benefit of using such characters in phylogeny. The mapping of ecological and behavioural features resulted in hypothesizing solitary behavior and life in forest as the conditions at the base of the mongoose clade.